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ABSTRACT 

After giving a pedagogical review of the chiral gauge approach to 2D gravity, with particular 
emphasis on the derivation of the gravitational Ward identities, we discuss in some detail the 
interpretation of matter correlation functions coupled to gravity in chiral gauge. We argue that 
in chiral gauge no explicit gravitational dressing factor, analogue to the Liouville exponential 
in conformal gauge, is necessary for left-right symmetric matter operators. In particular, we 
examine the gravitationally dressed four-point correlation function of products of left and right 
fermions. We solve the corresponding gravitational Ward identity exactly: in the presence of 
gravity this four-point function exhibits a logarithmic short-distance singularity, instead of the 
power-law singularity in the absence of gravity. This rather surprising effect is non-perturbative 
in the gravitational coupling and is a sign for logarithms in the gravitationally dressed operator 
product expansions. We also discuss some perturbative evidence that the chiral Gross-Neveu 
model may remain integrable when coupled to gravity. 

* unite propre du CNRS, associe a I'Ecole Normale Superieure et l'Universite Paris-Sud 
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1. Introduction 



What happens to a renormalizable two-dimensional field theory when it is coupled two 
gravity? In general we don't know. Of course, during the past years, tremendous progress 
has been made on a multitude of particular models, either through discrete matrix model 
techniques or in the continuum using the Liouville theory to describe gravity in the conformal 
gauge. Often, as is the case for the Ising model in a magnetic field, the coupling to gravity 
simplifies the theory alowing for an exact solution otherwise not available. Not withstanding 
these successes, the continuum methods are mainly restricted to conformal field theories. Also, 
we only have a very limited knowledge about correlation functions beyond the two- or three- 
point functions. Recently, by studying general continuum non-conformal field theories coupled 
to gravity in chiral gauge [1] it was shown that the one- loop /3-function gets affected by a 
universal gravitational factor. 

We think that the chiral gauge approach deserves further exploration. Here we report 
on some rather surprising results obtained in chiral gauge concerning a four-point function. 
This note is organized as follows: first, in section 2, we review the chiral gauge approach of 
Polyakov et al. [2-4] to 2D gravity, hopefully putting it into a pedagogical setting (see also 
ref. 5). In particular, we show how one obtains the gravitational action and the gravita- 
tional Ward identities. Then, in section 3, we discuss in some detail the interpretation and 
relevance of non-integrated correlation functions in chiral gauge. In particular, we point out 
why the non-integrated chiral gauge two-point functions give the gravitational scaling dimen- 
sions of the integrated conformal gauge two-point functions. We further argue that, in chiral 
gauge, for left-right symmetric matter operators Om no explicit gravitational dressing factor, 
analogue of the Liouville exponential of conformal gauge, is necessary, and that correlators 
like J d 2 x\ . . . J d 2 x n (Oyi(xi) ■ ■ ■ Om(^)) are well-defined quantities. Then, in section 4, we 
describe a sample computation of a four-point function using the Ward identities. It is first 
obtained as a perturbative solution of a partial differential equation. The perturbation series 
diverges but can easily be resummed. The resulting function, shown to be valid independent 
of perturbation theory, shows strong non-perturbative effects, namely logarithmic rather than 
power-type short distance singularities, invalidating the weak-coupling interpretation. Finally, 
in secttion 5, we discuss the relevance of this result to the gravitationally dressed operator 
product expansions, as well as the implications for the gravitationnal dressing of the inte- 
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grable chiral Gross-Neveu model. 



2. A review of 2D gravity in chiral gauge 

2.1. The matter action 

To start with, consider the action of a Majorana fermion \ coupled to gravity. Of course, 
one could consider a more general matter action as well, but let's be specific. 

Sm = ^J d 2 x(dete)x7 a e/«9 M x ■ (2.1) 

Here e aM is the zweibein and ej 1 its inverse. Our conventions are fairly standard* One finds 

S M = d 2 x [x-(-e + -d+ + e ++ d-)x- + x+(-e_+<9_ + e__<9+)X-] 

1 r / x / x i (2-2) 



V- (0+ - -±±<9_ j V- + ^ + (^cL - — d+ j *P + 

where we have rescaled ip+ = y/— e yx+ and ijj- = y/—e] X- (the local Lorentz phase can 

be chosen such that the square-roots are real). Note that although x+, X- are diffeomorphism 
scalars, ip + ,ip- bahave as half-differentials. 

If one now makes the conformal gauge choice e ++ = e = 0, e_| = e |_ = — so 

that (?-| = g |- = — e 2 ^ and g ++ = g = one obtains the well-known free-fermion action 

in conformal gauge: S'm = / d 2 x [ip-d+ip- + ip+d-tp+l. Here, however, we make a different 
gauge choice leading to the so-called chiral gauge: 

e+_e_ + = 1, e__ = 0, — = h ++ (2.3) 

e+- 

where the last equation is not a gauge choice but just the definition of h ++ . Then the fermion 



* Our conventions are: x ± — ± x 1 ), d± — -^{do ± di), x a y a = x + y + + X-y . The Minkowski 
metric is t] q = -1, r/n = 1 => rj+- = = -1, and g^ v = e ajl e a v , and e/e a v = 5» where, as 
usual, Lorentz indices (a, b, . . .) are raised and lowered with rj a b while is used for curved space indices 
(//, v, . . .). One defines \ = X + 7° = X7° an d 7°7 1 = 75 so that X7° = — x an d X7 1 — X75- Furthermore 
( 7 °)+ = -7°, ( 7 !)+ = 7 1 ,7 5 + - 75 and X ± = £(1±t&)x, so that X7 ± = -^xh°±J V ) = -V^Xt- Finall y 
note that for a 2 x 2-matrix M one has det M = M 00 Mu - M iM w = M ++ M - M+_M_+. 
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action becomes 



S M = J d 2 x hM0+ - h ++ d_)^ + ^J + d^ + ] . (2.4) 

Gravity is represented by the field h++ and only the left fermion ip- couples to gravity. It is 
straightforward to compute the metric tensor in this gauge: 

g+- = g-+ = g + ~ = g~ + = -i, g— = g ++ = o, g++ = -2h ++ , g — = 2h ++ (2.5) 

and —g = detg^ = 1. It is also easy to find that the only non- vanishing Christoffel symbols 



are 



r;_ = r: + = d-h++, r+ + = -d-h++, r; + = d+h++ + 2h++d-h ++ . (2.6) 

The stress tensor is defined in general as 

J__SS_ _L_ 8S 



and one obtains for the matter part in chiral gauge 



(2.8) 



It is straightforward to show that is classically conserved: V M T^ = 0. For v — + e.g., 
using (2.6) this is equivalent to the vanishing of 

0+T M + + cLTf + - (d+h ++ )T M _ - 2(d-h ++ )Tf_ - 2h ++ d-T M + , 

which in turn is shown using the equations of motion 

= h++d^ + ±(d-h++)i/>- , d-ij+ = . (2.9) 

Note also that the stress tensor is traceless: (T M ) M M = 0. 
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2.2. DlFFEOMORPHISMS 



Next, let us consider the effect of diffeomorphisms. Under a general infinitesimal diffeo- 
morphism — > x !± = + e ± (x + , x~) one has 

5e a ^(x) = e' (x) - e a ^(x) = e x d\e afl + e a \d^e x 

(2.10) 

8g^(x) = fi>(^) - 9^{x) = V^e v + = g up V ^t p + g^ v e p 

since V^^p = 0. Residual diffeomorphisms preserving chiral gauge must obey 5g = 5g + - = 

which implies (we write ej^ for residual diffeomorphisms) 

V^ = V_4 = & 8^ = 8-6+ = 0. (2.11) 

It is a straightforward excercise to check, using 5ip± = €^d\ip± ± ^(d+e^)tp±, 5h ++ = 
^d\h ++ + + 2/i++<9+e^, that the matter action (2.4) is invariant under residual dif- 
feomorphisms. We will also need to consider non-residual diffeomorphisms with general e~, 
but e + = 0. In this case one has to remember that h++ is defined as h++ = while 

g++ = —2e ++ e |_ = — 2(e_| e \-)h ++ . Although e |_ = 1 in chiral gauge, one has 

<5(e_| e |_) = cLe - (but still 5e = 0) so that Sg++ = — 2d-e~h ++ — 25h ++ . Thus 

Sh ++ = d+e~ + e-d-h ++ - h ++ d-e~ for e+ = . (2.12) 

Furthermore, for such diffeomorphisms 

5i[)+ = e~d-i[)+ 

1 _ for e + = (2.13) 

Sip- = e d-ip- + -(d-e 

so that one can show that the matter action (2.4) in chiral gauge is still invariant under these 
non- residual diffeomorphisms with e + = 0. 
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2.3. The gravitational action 



We want to compute the effective action T in chiral gauge, much along the lines of Friedan's 
Les Houches lectures for the conformal gauge [6]. The matter part is rather trivial since the 
original action is quadratic in the matter fields, but the gravitational part will turn out to be 
much less trivial. Start with the generating functional 

Z[V+,V-] = f^J ZWexp (-Sf(g) - S M (g,4>) + J d 2 x V^W-V- + (2-14) 

where Sff^g) — hq j d 2 x^/^g and VL&r is the volume of the diffeomorphism group. One has 

T>g = T>g++Dg T>g-\ = T>v + T>v~'Dh++dm x J where v are vector fields parametrizing the 

diffeomorphisms, m stands for the moduli, and the Jacobian J can be represented as usual by 
an integral over ghosts [3] 



% = J d 2 x [e++V-e- + C(V+e- + V_c+)] 



(2.15) 



so that 



z lv+,V-]=f dmVh ++ Z[ri + ,r]-,h ++ ,m] 
Z[ri + ,r]^,h ++ ,m] =e~ S « r (^J Vghe~ S ^ (^J V^e~ S ^ (2.16) 
xexp ^~ J (r ) -D+ 1 r)-+r) + dI 1 r) + )^ 

where D + = d + — h++d- and where we shifted as usual ijj- — > ip- + ^D^rj- and ip + — > 
i) + + \dZ X r] + . Let £ gh (#) = -log/ X>ghe _ ^ h and £ M (#) = -log/ Vi>e~ Sm . Then the 
generating functional of connected ^-correlation functions is* 

W[r] + , rj-, h++] = log Z[r] + , ??_, h++] 

= -sf - s gh (/ i++ ) - s M (^ ++ ) -\j {v-D+'v- + v+di 1 v+) ■ (2 ' 17) 



* Wc do not explicitly write the dependence on the moduli any longer. 
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To obtain the effective action T one introduces the classical fields ip± = |^ so that 

r[V>+,V>_,/i ++ ] = J d 2 x (t7_^_ + 77+^+) -W[r}+,ri-,h ++ ] 

= Sl Y + E gh (/i ++ ) + £ M (/i ++ ) + J d 2 x (iJ-D + ^ + ^+<M+) ( 2 - 18 ) 

= r [/i ++ ] + r exc it 

which we seperated into a (non-trivial) ground-state contribution To = + S g h + Sm and a 
trivial excitation part r exc it which is formally identical to the original matter action Sm ■ The 
effective action is to be used to compute correlation functions as 

$±(x 1 )..4±(x n )) = J dmVh ++ V^+V^ ^ ± (a;i)...^ ± (a; n )e- r ^+^-' ft ++] . (2.19) 

We have gone through the usual field theoretic formalism to show that the gravitational 
part To of T gets contributions from the ghost and matter sectors. Equivalently one could have 
argued a la David-Distler-Kawai [7] that in (2.16) the measures Dgh = V g gh. and Vijj = T) g ip 
are complicated but factorize into e~ Sgh r>ogh and e~ SM r , oV ; where £>ogh and T>Qip now are 
trivial (flat) measures. Then one would directly obtain (2.19) (except for the replacement 

V>± -> V±)- 

Now, if the Faddeev-Popov procedure of factorizing the volume of the diffeomorphism 
group is to make sense, then T must be invariant under general diffeomorphisms, not only 
residual ones.^ To check the invariance of T properly one should have imposed, instead of 

(2.3), conditions like e_| e |_ = a, §=■ = (3 and the definition = h ++ , and keep a and 

(3 as classical non-dynamical background fields. Then under a diffeomorphism (2.10) one has 
e.g. 5(3 |^ _ Q = d-e + . The classical action 5m would then read 

S M = J d 2 x i^(d+ - h ++ d_)^ + ^ + («9_ - [3—d + )^+] (2.20) 

which is obviously invariant under diffeomorphisms, since it is just the original invariant matter 
action written in a particular form. Actually, for our purpose of deriving the gravitational 

f One may view this requirement as a condition on the counterterms. 
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action To is is enough to consider diffeomorphisms with e + = and arbitrary. In this case 
there is no need to introduce the a and /3 into Sm and one can directly work with Sm as given 
by (2.4). As already noted earlier, (2.4) is invariant under diffeomorphisms with e + = 0. Since 
Sm is formally identical with r exc i t the same is true for the latter. Thus imposing invariance 
of T translates into imposing invariance of To [g^] : 

= 6T = J d 2 x5g^^ = J d 2 xv^2VWs r = _ 2 J ^ x ^- ge ^ T ^ (2 . 2 1) 
where 

T grav _ \ icy 99 \ 

Equation (2.21) implies the conservation of the gravitational stress-energy tensor, in particular 
with e+ = 0: 

= v^T grav = -V+T grav + 2/i ++ V_T grav - V_T grav 

1 (2-23) 
= V+T grav + h++ V_T grav + -V- (T grav )^ 

or 

[d + - h ++ d_ - 2(d_h ++ )] T grav = iv_(T gra T M . (2.24) 

Now one uses the fact that (T grav ) M M is a gravitational scalar of dimension two, and hence 
must be proportional to the curvature scalar: 

(T grav )% = • (2.25) 

Combining eqs (2.24) and (2.25) gives 

[d+ - h ++ d_ - 2(d-h ++ )} T grav = ^d-R (2.26) 

which is the central equation for determining the gravitational action Tq. It is straightforward 
to compute R in chiral gauge: 

R = 2<9 3 . (2.27) 
Furthermore, since ^/^g — 1 one nas T grav = \ and (2.26) turns into a functional differ- 
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ential equation for Tq: 



[d+ - h ++ d- - 2(d-h ++ )} ^±±1 = ^<9 3 h++ . (2.28) 

Its solution is 

= »/d», ^-h ++ ) g _ {g+ ! h++a _f -h ++ + „ /d»* . (2.29) 

Indeed, observing that (d-d + — d-h ++ d )~ 1 is a symmetric pseudodifferential operator one 
immediately finds 

Applying then d + — h++d- to the r.h.s. of this equation yields ^d^_h ++ + 2<9_/i ++ T grav , 
showing that (2.28) is satisfied. Since d-(d+ — h++d-) is — \ times the Laplacian on a scalar, 
one finds of course that To can be written as 



r ° = ~48j dx ^~ aR V2 R + VJ dx V=9 (2-31) 

which is the famous covariant form written by Polyakov [3], and which, in conformal gauge, 
reproduces the Liouville action. 

2.4. Gravitational Ward identities 

Note that equation (2.28) can be rewritten, using (2.12), as 



/ 



d 2 x5h ++ ^^ = / d 2 xe-d 3 _h ++ . (2.32) 
oh-\ |_ 12 



This means that To, if considered as a functional of h++ only, is not (completely) invariant 
under diffeomorphisms (even with e + = 0) but has the anomaly as given by the r.h.s. of 
this equation. Equation (2.21) however, expresses the diffeomorphism invariance of To as a 

covariant functional of g ++ and (and g ), where the chiral gauge is only imposed after 

the variation is performed. 
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It is now straightforward to derive Ward identities for correlation functions: 



{ip(x 1 )...ip(x n )) = J Vh++e- r ° [h++] J ViPiP( Xl )...ip(x n )e^ cit[h++ ^ ] . (2.33) 

Here ip need not be the previous fermion fields, but could be more general matter fields. 
Changing variables of integration to h ++ = h ++ + Sh ++ and ip — ip + Sip such that Sh ++ and 
Sip correspond to the diffeomorphisms (2.12) and (2.14), hence 5r exc i t = 0, one obtains, using 
(2.32), the following Ward identity 

■ ■ ■ ■ • • V>M> + —J d 2 ze-(z)(& i _h ++ (z)iP(x 1 ) . . . iP(x n )) = . (2.34) 



Next, following ref. 2, we will turn this Ward identity into a partial differential equation 
for (ip(x{) . . .ip(x n )). To do so one has to eliminate h++ from the correlation function. Here, 
we will again concentrate on the above example of fermions. The classical equation of motion 

d + ip_ = h ++ d-ip- + A(d-h ++ )ip_ (2.35) 

with A = \ carries over to the quantum case but for two differences: first, the product of 
fields at the same point needs to be regularized by some normal-ordering prescription which 
may modify the weight A. Second, when inserted into the functional integral, the equations 
of motion remain true up to contact terms (e.g. ( §^\ <P(y)) — S(x — y)). Covariance of the 
quantum equation of motion (2.35) requires that the variation of ip- under diffeomorphisms 
also gets modified: eq. (2.13) gets replaced by 

Sip- = e~d-ip- + A(d-e~)ip- (2.36) 

with the same A as in (2.35). The Ward identity then becomes 

.(xi) . ..ip-(x n )) 



\9 d 3 

+ ^ Q(p^ ( h ++( z w-M ■ ^-M = o . 



(2.37) 
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Using the identities 



<9 3 [z - x i f 



d{z~f z+ 



d 3 2(z~ - x~) 
d(z') 3 z + —xf 

eq. (2.37) can be integrated as 



ATxi5^\z - Xi) , 
d 

4m —5 {2) (z - Xi) , 



(2.38) 



E 

i=i 



d 



z T — x A 



dxi 



2A- 



Z T — Xi 



iirX 9 



(h ++ (z)il)-(xi) . ..ip-(x n )) . 



(2.39) 



Now one uses the quantum equations of motion (2.35) as: 



d 
dz+ 



(l/j_(z)l/j-(x2) . ..ifj-(Xn)) 



= ((h ++ (z)d^(z) + Ad-h ++ (z)lP-(z)) V-(X 2 ) • • • V-M) 

= ^{h ++ (z)^(z)^(x 2 ) . . . + (A - l){(d-h ++ (z))lP-(z)lP-(x 2 ) ■ ■ ■ 1>- (X n )) 



(2.40) 



(up to contact terms). Using (2.39) and its derivative w.r.t. z one finally arrives at 
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d 

dz+ 



i=2 



(z — Xi ) 2 d d A z — Xi 
+ 2 A 



z+ - x\ dz dx i 



d d 



2A< 



: " - xt V dx i dz 

<lM*i). 



z + — xt 



(2.41) 



where we set 



7 = 



(2.42) 



We quote without proof [3] that 7 is related to the total central charge c of the matter coupled 
to gravity (e.g. c = | for a Majorana fermion) by the relation 



7 = !^( C ~ 13 ~ V^-lXc-25)) . 



(2.43) 



Equation (2.41) was first written in ref. 1. 
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As an example for the use of the Ward identities, consider the two-point function of ip-. 
Perturbation theory* suggests the ansatz 

// f\i f \\ 1 [( x ~ -y~)( x+ -y + )T 25 (OAA , 

W-(xW-{y)) ~ — . 1+2S f + = = = • ( 2 - 44 ) 

Gravity only contributes a left-right symmetric factor [(x~ — y~)(x + — y + )] . Furthermore, 
by the usual arguments, \ + 5 must coincide with the anomalous dimension A of eq. (2.36). 
Inserting this ansatz into (2.41) for n = 2 gives an algebraic equation for A: 

. A(A-l) 
A - A = — (2.45) 

7 

where Ao = \. This is the well-known KPZ-equation [4] expressing the anomalous dimension 
A in the presence of gravity in terms of the dimension Aq without gravity. 



3. Interpretation and relevance of 
non-integrated correlation functions 

What do correlation functions like (2.44) mean? Since one has integrated over the metrics, 
i.e. over h++, what is the distance between x and y? Clearly, these are non-trivial questions. 
Let us compare with what one does in conformal gauge. In conformal gauge, one usually 
computes integrated correlation functions at fixed area A, like 



( / d 2 xe a ^O(x) [ d 2 ye a ^O(y)} (3.1) 

J J fixedA 



where O is a left-right symmetric matter field of conformal dimensions (Ao,Ao) (e.g. the 
product of our fermion fields ip- and ip + with Ao = \) and the Liouville field. The constant a 
is chosen such that the conformal dimension of e a ^ is (1— Ao, 1— Ao), so that the total integrand 
has conformal dimensions (1, 1), and the integral is invariant under conformal transformations. 
The area of the surface can be fixed by adjusting the zero-mode of the Liouville field <fi. 



* One can do a simple Feynman diagram expansion of the two-point function using the vertices and 
propagators derived from T. This is a perturbation scries in ~ i. 
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Correlation functions like (3.1) are conformal scalars, i.e. are invariant under the residual 
diffeomorphisms of conformal gauge and have a well-defined meaning. It has been shown [7] 
that (3.1) scales with the area as ^4 2 ~ 2A where the gravitational scaling dimension A is given 
by the KPZ formula (2.45). Hence this A coincides with the A-exponent characterizing the 
non-integrated two-point function in chiral gauge: 

(0(x)0(y)) ~ (x __ y _ )2A( ^_ y+)2A _ 2Ao >< (x+ _ 1 y+) 2A ( 3 - 2 ) 

where e.g. in the case of the fermions the first factor comes from ■?/>_ and the second from 
■0+. One sees that although such non-integrated correlation functions in chiral gauge are not 
invariant, their singularity structure (exponent A) nevertheless has an invariant meaning. Let 
us try to understand why A as given by (3.2) should coincide with the gravitational scaling 
dimension of (3.1). In chiral gauge, since ^/~# = 1> the area of a surface is completely 
independent of the metric h++. Whereas in conformal gauge one could choose the range of 
the coordinates x + ,x~ to be fixed, in chiral gauge their range is relevant to the geometry 
(and is part of the moduli of the surface). Integrating (3.2) in x and y over the surface 
then gives A 2 ~ 2A , possibly up to an A-independent constant. Thus the A characterizing the 
power-law behaviour of the non-integrated two-point function directly gives the gravitational 
scaling dimension without further dressing by some field f(h++) that would be the chiral 
gauge analogue of the e a< ^-dressing. To understand why no such extra dressing is required in 
chiral gauge, let's go back to the example of the fermion fields, i.e. O = ip+ip-. Under an 
e~-diffeomorphism we had 

otp+ = e £/-Vm- 
and similarly one finds for an e + -diffeomorphism 

6if>_ = e+<9 + ^_ + (A - -)(d+e+)^ 

2 (3.4) 
= e+d + ^J + + -(d+e+)ij + 

which combines into 

= e X d x ^ + A(<9 A e A )V- - l(d+e+)^ 

(3.5) 

^ + = e A ^ + + i(<9 +e +)^ + . 
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For residual diffeomorphisms (preserving chiral gauge) one has <9a 6 r = and it follows 

<5rV± = 4d x if>± ± \(d+e + )^ ± (3.6) 

and as a consequence 

M^MO = e^AftMO = d x (4M-) (3.7) 

which is a total derivative. Hence J d 2 iO(i) = J d 2 x ^jj + (x)^jj-(x) is invariant under residual 
diffeomorphisms, which is the analogue of the (1,1) condition in conformal gauge. We see that 
no extra gravitational factor f(h++) is needed to achieve this, explaining why the exponent 
A = | + 5 of (2.44) directly gives the gravitational scaling dimension. 

This result is of course not restricted to fermions. From (3.3)-(3.7) we see that if under a 
general diffeomorphism 

SO = e x d x O + A{d x e x )0 (3.8) 
then under a residual diffeomorphism 

SrO = eid x O = <9 A (e R 0) =s> 5 R Jd 2 xO = 0. (3.9) 

Thus integrated n-point functions f d 2 xi . . . d 2 x n (0(xi) . . . 0(x n )) are perfectly well-defined* 
objects. Of course, it is just as meaningful to consider also the non-integrated n-point 
functions which are scalar densities w.r.t. residual diffeomorphisms. The invariance of 
the integrated n-point functions under residual diffeomorphisms should be expressible as 
BRST-invariance.^ To our knowledge, the BRST operator for the chiral gauge has not yet 
been constructed, nor its cohomoly been investigated, but it is clear that correlators like 
(J d 2 xi?p + (xi)ijj-(xi) . . . J d 2 x n ip+(x n )4>-(x n )) should turn out to be BRST-invariant. 



* Let us insist that well-defined means invariance (or covariance) under residual diffeomorphisms. Since 
one works with a fixed gauge this is as much as one can demand. Of course, one should be able to compare 
to another, say conformal gauge, but this is beyond the scope of this note. 

f Recall that in conformal gauge, BRST-invariance is essentially the statement that the integrand is a 
(1, 1) field w.r.t. conformal transformations, i.e. residual diffeomorphisms. 
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4. The fermion four-point function 



As a non-trivial example of how to use the Ward-identity (2.41) we now compute the 
fermion four-point function. As discussed in the next section, this is also of some relevance to 
the gravitational dressing of the chiral Gross-Neveu model. Therefore we also add some colour 
indices i,j — l,...,N, and take as the matter action N copies of (2.4). We want to compute, 
for i 7^ j 

G 4 (w,x,y,z) = (^L(«;)^(ti;)^L(x)^(x)^(y)^(y)^i(«)^.(«)> 

1 



(w + — x + )(y + — z + ) 
In analogy with the two-point function we use the following ansatz 



where the anharmonic ratio t = t~ is given as usual by 

(w~ — y~)(x~ — z~) , t . 

t = t = _ _ _ rr 4.3 

{w - z ){x -y ) 

and similarly for t = t + . Inserting the ansatz (4.2) into the Ward identity (2.41) with n = 4 
leads after some algebra to* 



7idt + ^—l(t-t)d t td t + (1 - AA)td t + 2A 2 ^| 



f(t,t)=0, (4.4) 



Note that A = 1 reproduces the equation derived in ref. 1 for the four-current correlation 
function. The ansatz (4.2) is justified by the fact that we obtain an equation for f(t,t) 
involving only t, i, dt, df and not x, y, z or w explicitly. 

The partial differential equation (4.4) has many solutions. To pick out the physical ones 
we have to compare with perturbation theory in ^. Looking at the effective action T (cf. 
(2.18), (2.29)) it is clear upon rescaling h ++ = ^h ++ that each interaction involving a h ++ 



| Performing the change of variables carefully also leads to contact terms involving 8^(x — y) or 5^ (w—z). 
We will assume that x ^ y and w ^ z, so that we can drop these contact terms. In particular, below we 
will consider t —* 1, i.e. nioori/^z which is perfectly compatible with this assumption. 
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is accompanied by at least one factor of Doing a standard perturbative Feynman diagram 
expansion, we obtain 

1 / 4- 1 2A 2 / 4- 1 

/ = 1 - --— \ogti+ 0(lh 2 ) = 1 -—- \ogti + 0(l/ 7 2 ) (4.5) 

Z'-f t — 1 7 t — 1 

since A = ± + 0(1/7). 

In principle, the partial differential equation (4.4) could be solved order by order in a 
perturbation series in ^. In practice, this leads to very complicated poly-log integrals already 
at low orders. The main difficulty is the factor 1/(1 —t). If one considers the vicinity of t — 1, 
this difficulty disappears, and (4.4) becomes 



-ytdt + {t- l)d t td t + (1 - 4A)<9 t + 1 1 " 



t 



/i(M~)=0 (4.6) 



where the subscript 1 on / is to remind us that f\ ~ / only in the vicinity of t — 1. This equa- 
tion (4.6) can be solved exactly. Writing the solution as a perturbation series in ^ (matching 
to (4.5)) gives 

h{t,t) = h{9)-2^—-\ 9 , 9-~ ,._ ^ (4-7) 

n=0 ' ' V y 

where a n = a(a + l)(a + 2) . . . (a + n — 1). At each order in ^ one can of course replace logt 
by logtt in g, as suggested by (4.5). The series (4.7) has zero radius of convergence, but its 
Borel transform can be recognized as the hypergeometric function 



B[fi](u) = £ I^fC u n = F(2A, 2A, 1; u) . (4.8) 

71/. 71/. 



n=0 

The inverse Borel transform 



00 

h(g) = J dve- v B[h]{uv) (4.9) 


gives the resummed function f\(g) in terms of a Whittaker function. Alternatively, one can 
directly observe that (4.7) coincides up to an overall factor with the asymptotic expansion of 



Then /1 is no longer an exact solution of (4.6) but an exact solution of another equation, differing from 
(4.6) only by higher order terms in (i — 1). 
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a Whittaker function. Hence 



/ ,\2A 1 / A 2A-l/2 1 

= ("J *(2A, 1; --) = (--J e- 1 ^ W 1/2 _ 2A>0 (-_) . (4.10) 

Here W is the Whittaker function and ^ is a solution to the degenerate hypergeometric 
equation [8]. Indeed, although the equation (4.6) has many solutions, if one uses an ansatz 
with fi only depending on t and t through g, then equation (4.6) becomes 

{<? 2 ^ + [(1 + 4A)s - 1] A + 4A 2 } h{g) = . (4.11) 

Setting fi(g) = (—l/g) 2A u(—l/g) one sees that u(x) satisfies the degenerate hypergeometric 
equation 

xu"(x) + (b - x)u'(x) - au(x) = (4.12) 

with b = 1 and a = 2 A. Perturbation theory has told us which of the two independent 
solutions to choose, namely u(x) = \I/(2A, l;x). Let us insist that we just showed that (4.10) 
is a solution to the differential equation (4.6), independent of perturbation theory in ^. 

Having the perfectly non-perturbative expression (4.10) for fi(g), we can now investigate 
its behaviour for large g, which is just the series expansion of \I/(2A, 1; x) for small x: 



k 



-I) (4,3, 



where tp(x) = T'(x)/T(x). Recall that g = — \ogi/[y(t — 1)], hence large g means t — > 1 (for 
fixed t), so this is the limit where / ~ /i- Although (4.13) is the exact asymptotic for /i, it 
gives only the leading order for / (i.e. we can only trust the k = term): 



f(t,t) 



2A 



7(*-l)Y 1 



logt J T(2A) 



^(1) - log ( ^-h ) +0 ((t-l),(t-l) log(f - 1)) 



(4.14) 

What does this mean for the fermion four-point function (4.1)? Since 

t _ fa- -'-Ihr-O (415) 

(u; -z )(x -y ) 

one has t — > 1 if either w~ — > or y~ — > z~ , i.e. when two fermion operators of the same 
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colour approach each other. Inserting (4.14) into (4.1) and (4.2) then gives 



2A 2 ^ 2 ^ 

G A (w,x,y,z) ~ Y(2A) ~ z ~^ x ~ ~ y ~^~ t lo S^ 2A [i w+ ~ x+ )(v + ~ z+ )] ~ 



x 



, log A , (w — x )(y — z ) 
^(l)+log -M - log 1 



7 J (w -z )(x -y ) 



(4.16) 

It is important to realize that we work in Minkowski space so that we can take t — > 1, 
keeping t ^ 1 fixed. Rather surprisingly, the four-point function (4.16) no longer contains 
the perturbative singularity ~ (w~ — x~)~ 2A (y~ — z~)~ 2A , but resumming the series has 
transformed it into a logarithmic singularity, plus a non-singular part!* 

Mathematically, the origin of the logarithm can be traced to the degenerate hypergeomet- 
ric equation (4.12) satisfied by u(x) = x~ 2A fi(—l/x). For generic parameter b it has two 
independent solutions [8] $(a, b;x) and x 1 ~ b ^(a + 1 — b, 2 — b; x). Obviously, for b — > 1 the 
second solution generates log a; $(a, b; x), among others. This is a well-known phenomenon in 
the theory of ordinary linear differential equations. 

Physically however, it was quite unexpected that turning on gravity (^ ^ 0), even infinites- 
imally weakly, completely changes the singularity structure: this is a truely non-perturbative 
phenomenon, due to the divergence of the perturbative series in K 



5. Conclusions and Outlook 

What do we learn from all these computations of the four-point function? One obvious 
lesson is - unlike the situation of the two-point function - that we cannot trust the weak- 
coupling gravitational perturbation theory in ^ ~ What is the meaning of the (non- 
perturbative) logarithmic singularity of the four-point function? In section 3, we have argued 
that integrating correlation functions like (4.16), computed in chiral gauge, leads to well- 
defined objects, invariant under residual diffeomorphisms. Obviously, we cannot integrate 
our result (4.16) since it is valid only in the vicinity t ~ 1. However, since integrating is 

★ Naively it looks as if (4.16) now contains a new singularity ~ (w~ — z~)~ 2A (x~ — y~)~ 2A as w~ — > z~ 
or x~ — > y~ . However, this means t— 1 — > oo which is clearly outside the domain of validity of eq. (4.16). 
Let us also insist that the contact terms ~ 6^ 2 \x — y), 8^(x — y) we dropped above precisely correspond 
to t — 1 — > oo and are completely irrelevant to the behaviour as t — > 1. 
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a "trivial procedure" it certainly makes sense to study the properties of the non-integrated 
correlation functions as well. The main question that arises is whether the correlator (4.16) 
tells us something about the gravitationally dressed operator product expansion. Equation 
(4.16) should express the OPE of (j>(w) = ip-(w)ip + (w) with (f>(x) = ip^(x)ip + (x) as w~ — > x~ . 
It looks like 



(w)(f>(x) ~ 0(x) + log(u> — x )0(x) . 



(5.1) 



This is a particular example of a more general OPE with logarithmic short- distance behaviour: 



(w)(f)(x) ~ - x) An ~ 2A4, O n + . . . + log(tU - x)O n + ... 



(5.2) 



In conformal gauge, such logarithms have been noticed before in ref. 9 where the WZW 
model based on the supergroup Gl(l, 1) was discussed. Later a more systematic discussion 
was given in ref. 10, where the appearance of logarithms in conformal blocks in the c = —2 
and other non-unitary models was studied. There it has been argued that the emergence of 
logarithms in correlation functions is due to new so-called logarithmic operators, whose OPEs 
display logarithmic short-distance singularities. These new logarithmic operators have confor- 
mal dimensions degenerate with those of the usual primary operators, and it is this degeneracy 
that is at the origin of the logarithms (cf. our discussion of the degenerate hypergeometric 
equation in the previous section). As a result one can no longer completely diagonalize the 
Virasoro operator Lq, and the new operators, together with the standard ones form the basis 
of the Jordan-cell for Lq. In the case of two operators with degenerate conformal dimensions 
A n the operator product expansion precisely takes the form (5.2), while the OPE with the 
conformal stress-energy tensor is 

T(z)O n (0) ~ ^O n (0) + \o n (0) + -dO n (0) , (5.3) 

in particular 

L \O n ) = A n \O n ) , L \O n } = A n \O n ) + \O n ) . (5.4) 

This makes it possible to have logarithmic terms in the correlation functions without spoiling 
the conformal invariance. 
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Finally we would like to comment on the relevance of our present results for the gravi- 
tational dressing of a two-dimensional integrable but not conformally invariant field theory, 
namely the chiral Gross- Neveu model [11]. As is well-known, its action is given by the massless 
free-fermion action (2.1), where the fermions are N component fields, and an interaction term 
between two left and two right fermions ~ J d 2 £ In chiral gauge, only the left 

fermions ^_ interact with gravity, cf. eq. (2.4). Without gravity, it is known that this model is 
completely integrable [12] and exhibits dynamical mass generation [11]. Does the integrability 
remain once the model is coupled to gravity? A necessary condition of integrability is that the 
S-matrix for the scattering of the physical particles (here the massive fermions) is factorizable 
and elastic. This actually is a consequence of the factorizability and elasticity of the S'-matrix 
for the pseudoparticles (here the original massless fermions). 

As the simplest check, we have investigated whether the two-pseudoparticle S'-matrix re- 
mains elastic in the presence of gravity. Here again we face the issue of the interpretation of 
the S-matrix elements in the presence of gravity: the S-matrix, e.g. for the scattering of two 
left fermions, is obtained from the four-point function (4.2) by removing the external propa- 
gators and setting the external momenta on-shell (p + = p' + = q + = q' + = where p,p' and 
q,q' are the initial and final momenta). According to our discussion of section 3, this does 
not seem to lead to a well-defined quantity. However, bearing in mind the "experimental" 
situation for measuring S-matrix elements, even in the presence of gravity, we expect that 
the S-matrix should be well-defined at least within a gravitational weak-coupling expansion, 
provided the latter makes sense. Now for the chiral Gross-Neveu model, there is no scattering 
of two right pseudoparticles, while the left-right scattering is always elastic, as can be seen 
simply by combining momentum conservation and the on-shell condition. It remains to con- 
sider the scattering of two left fermions which interact due to their coupling to gravity, cf. 
eq. (4.2). To first order in ^, this S-matrix element vanishes. Indeed, from (4.5) one finds 
upon Fourier transforming that is is given by ^(p+p')-(q + q')- ^z^tw which vanishes on-shell 
(p+ = P+ = 0)- We have verified that this remains true at the next order in ^, including "two 
graviton exchanges". If this remains true at all orders in ^ and even non-perturbatively, one 
would have complete elasticity of the two-pseudoparticle S-matrix, and this would certainly 
be a well-defined and gauge-invariant statement. One could then go on and speculate that 
all S-matrix elements for the pseudoparticle scattering remain elastic and factorizable in the 
presence of gravity and that the same is true for the physical (massive) fermions, in other 
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words that the integrability of the Gross-Neveu model survives coupling to gravity. However, 
there is still a long way to go. 
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